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Abstract 

Kuchaf showed that the quantum dynamics of (1 polarization) cylindrical wave solutions to vacuum 
general relativity is determined by that of a free axially-symmetric scalar field along arbitrary 
axially-symmetric foliations of a fixed flat 2+1 dimensional spacetime. We investigate if such a 
dynamics can be defined unitarily within the standard Fock space quantization of the scalar field. 

Evolution between two arbitrary slices of an arbitrary foliation of the flat spacetime can be built 
out of a restricted class of evolutions (and their inverses). The restricted evolution is from an initial 
flat slice to an arbitrary (in general, curved) slice of the flat spacetime and can be decomposed 
into (i) 'time' evolution in which the spatial Minkowskian coordinates serve as spatial coordinates 
on the initial and the final slice, followed by (ii) the action of a spatial diffeomorphism of the 
final slice on the data obtained from (i). We show that although the functional evolution of (i) is 
unitarily implemented in the quantum theory, generic spatial diffeomorphisms of (ii) are not. Our 
results imply that a Tomanaga-Schwinger type functional evolution of quantum cylindrical waves 
is not a viable concept even though, remarkably, the more limited notion of functional evolution 
in Kuchaf 's 'half parametrized formalism' is well-defined. 
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I. INTRODUCTION 



Given the difficulties of constructing a quantum theory of gravity, it is of use to build 
intuition from the study of simpler toy models. While cosmological mini-superspaces 0, 
have yielded valuable insights, their finite-dimensional nature precludes the occurrence of 
field-theoretic aspects of quantum gravity. In this regard, vacuum cylindrically symmetric 
gravitational fields with 2 hyper surf ace-orthogonal, commuting Killing vectors (one gener- 
ating translations along 'Z' and the other, rotations along '$') constitute a useful infinite- 
dimensional midi-superspace. Quantization of this cylindrical wave midi-superspace was 
initiated by Kuchaf and further studied by a number of authors 0, S S S H|- While 
the latter works deal with quantization after fixing a gauge, our focus here is on aspects of 
Dirac quantization of cylindrical waves. 

Every cylindrical wave solution is determined by a corresponding axially-symmetric so- 
lution to the free scalar wave equation on a fixed 3 dimensional flat spacetime. Specifically, 
the cylindrical wave line element is 

ds 2 = e^i-idT) 2 + (dR) 2 ) + R 2 e~^(d^) 2 + e^(dZ) 2 (1.1) 

and that of the flat spacetime is 

ds 2 = -(dT) 2 + (dR) 2 + R 2 (d$) 2 . (1.2) 

ip{R,T) is the free scalar field and ^y(R, T) is its energy in a box of size R at time T. 
Kuchaf defined a canonical transformation from the Arnowitt-Deser-Misner (ADM) phase 
space of cylindrical waves to that of an axially-symmetric parametrized field theory on 2+1 
dimensional flat spacetime and studied the formal Dirac quantization of the latter system 

1- 

Parametrized field theory (PFT) is free field theory on flat spacetime in a diffeomorphism 
invariant disguise Q. It describes field evolution on arbitrary (and in general, curved) 
foliations of the flat spacetime instead of only the usual flat foliations, by treating the 
'embedding variables' which describe the foliation as dynamical variables to be varied in the 
action in addition to the scalar field. In the present context the coordinates X a := (T, R) are 
parametrized by a new set of arbitrary coordinates x a = (t, r) such that for fixed t, X a (t, r) 
define an axisymmetric spacelike slice in Minkowski spacetime with radial coordinate r. 
General covariance of PFT ensues from the arbitrary choice of x a . In the Hamiltonian 
description, the embedding variables l X al are canonical coordinates and general covariance 
implies that Hamiltonian evolution from one slice of an axisymmetric foliation to another is 
generated by constraints. 

Indeed, under Kuchaf 's canonical transformation 1 , appropriate combinations of the 
Hamiltonian and diffeomorphism constraints of the midi-superspace take the form of the 
constraints appropriate to PFT, namely, 

C a {x) :=P a (x) +h a [ip,n,X a ](x) « , (1.3) 

where P a and it are the momenta conjugate to X a and ip, an d h a is related to the stress- 
energy of the scalar field 2 . The gauge fixing conditions X a = x a 3 map the classical 



1 While R is directly related to the radial components of the metric Kuchaf 's canonical transformation 
identifies T with the spatial integral of one of the extrinsic curvature components of the cylindrical wave 
metric (|l.lf) on a t— constant slice. 

2 As is standard in constraint dynamics, w denotes weak equality. 
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theory directly onto the standard description of a free scalar field in Minkowskian co- 
ordinates (R/T). This description admits the standard Heisenberg picture based Fock 
quantization [3J |4| . 

In the Dirac quantization, a formal operator version of C a acting on a physical state 
of the theory is given by 

Eq. (|1.4j) takes the form of a functional Schrodinger equation which represents infinitesimal 
evolution of the Schrodinger state from one Cauchy slice to another. In this work 
we investigate if the formal equation (J1.4)) is well defined in the context of the standard 
Heisenberg picture based Fock space quantization of the free scalar field on flat spacetime. 
More precisely, in analogy with the relation between the Heisenberg and Schrodinger pictures 
in quantum mechanics, we investigate if there exists a unitary transformation on the standard 
Fock space which implements evolution between any two axisymmetric Cauchy slices of the 
flat spacetime and whose infinitesimal version yields equation (jl.4j) . 

What do we mean by evolution between two slices? Since the classical theory is that of a 
free scalar field we explicitly (and uniquely) know the evolved data (tp, it) on the slice t = t\ 
in terms of initial data at t — to- Since the equations for scalar field evolution are linear, 
we also know how the corresponding quantum operators are related. This is what we mean 
by evolution. Since this notion of evolution is independent of the choice of interpolating 
foliation between the two slices, our subsequent considerations are phrased solely in terms 
of data on the two slices. 

We enquire if the operators at the initial and final slices are unitarily related. If they are, 
then in analogy to the definition of the Schrodinger picture from the Heisenberg picture in 
usual quantum mechanics, we can define the Schrodinger state at any time t as the unitary 
image of the state at time to 4 5 . Techniques appropriate to our investigation here have been 
developed in (IoLlll| for the case of PFT on (n+l)-dimensional flat spacetime in the absence 
of axisymmetry. That work arose as an effort to rigorously implement the considerations 
of Kuchaf in j§, The strategy pursued in [HI 11 is as follows. Classical scalar field 



evolution from an initial flat slice to any final slice is a linear canonical transformation on 
the scalar field phase space. Instead of phrasing this canonical transformation in terms 
of the scalar field and its momentum, one can instead use appropriate linear functionals 
of these fields. The linear functionals used are such that when they are evaluated on the 
initial slice, they reduce to the familiar mode coefficients (i.e. the classical correspondents 
of the standard annihilation and creation operators). It then follows (from the fact that 
evolution is a canonical transformation) that these functionals evaluated on the final slice 
are related by a Bogolubov transformation to their initial values. The criterion for unitary 
implementability of a Bogolubov transformation is that the '/?' matrix be Hilbert-Schmidt 
(i.e. the squares of the absolute values of its diagonal elements should be summable) [3, 1^ • 



3 Technically, this is not a gauge fixing but a 1-parameter family of gauge fixings (one for each value of t) 
known as 'deparameterization'. 

4 For a beautiful exposition of the Schrodinger and Heisenberg pictures in a canonical treatment of PFT 
on (n+l)-dimensional flat spacetime (in the absence of axisymmetry), we refer the reader to 0,0]. 

5 Note that if we restrict attention to only flat slices, the standard Hamiltonian is the generator of the 
desired unitary transformation. The non-triviality enters solely due to the possibility of evolution along 
directions which are not isometries of the spacetime. 
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The results of |10L lllj are as follows. In 1+1 dimensions operator evolution is unitary, 
the functional Schrodinger picture exists as the unitary image of the standard Heisenberg- 
picture-based Fock quantization and the infinitesimal version of the unitary transformation 
on Heisenberg states implies that the Schrodinger states satisfy an equation of the form (|1.4|) 
(albeit, with an extra 'anomaly potential' term In n+ 1 dimensions (n > 1), the (3 



coefficients for evolution along generic curved-foliations are not Hilbert-Schmidt. Hence the 
functional Schrodinger picture does not exist (at least in the form usually envisioned, as the 
unitary image of the Heisenberg picture) and the functional Schrodinger equation cannot be 
given any obvious meaning. The case of interest here, namely 2+1 dimensions with fewer 
degrees of freedom (only the axisymmetric ones), is a tantalizing intermediate case and as 
we shall demonstrate, we get a suitably 'in between' answer. 

Here, our interest is in evolution between two axisymmetric, but otherwise arbitrary, 
Cauchy slices Ej and Ey in the flat spacetime (|1.2j) . Note that the restriction of the space- 
time coordinate R to any axisymmetric spatial slice defines a natural radial coordinate on 
the slice. We refer to this natural coordinate by the same symbol R. In order to accommo- 
date arbitrary reparameterizations {R,T) —>■ (r,t), it is necessary to admit arbitrary radial 
coordinates Tf(R) on Ej and r^R) on Ej 6 . Despite this, it suffices for our purposes, to 
consider evolution from a fixed initial slice Eo, chosen to be the T = slice coordinatized by 
the natural radial coordinate R, to an arbitrary final slice E, coordinatized by an arbitrary 
radial coordinate r = r(R). This entails no loss of generality since, if in the quantum theory 
all evolution is unitary, quantum evolution from any Ej to any Ej may be constructed as 
the (inverse) unitary evolution from Ej to E followed by unitary evolution from E to Ej. 

Now, consider the restricted evolution from E to E such that the radial coordinate r on 
E is chosen to be the natural radial coordinate R i.e. r is identified with the restriction of 
the spacetime coordinate R to E. We shall refer to such evolution as the 'half parametrized' 
evolution pj. Next, we reconsider the general case of interest, namely that of evolution from 
So (with coordinate R) to E with coordinate r = r(R). The uniqueness of evolution from 
initial data implies that {ip{r), 7r(r)) on E can be obtained by first evolving the initial data 
from Eo to E by a 'half parametrized' evolution and then subjecting the result to the spatial 
diffeomorphism defined by r(R). 

In view of this we proceed as follows. In section II we construct the Bogolubov trans- 
formation corresponding to the most general evolution from E to E. In section III, we 
particularize the Bogolubov transformation of section II to the case of 'half parametrized' 
evolution defined above and show that the (3 matrix is Hilbert-Schmidt. This implies that 
the finite functional evolution corresponding to Kuchaf's 'half parameterized' formalism [3] 
is well defined. In section IV we show that the action of a generic spatial diffeomorphism on 
data on any slice E is not unitarily implemented in quantum theory. This implies that the 
most general Tomonaga- Schwinger type functional evolution of quantum cylindrical waves 
is not a viable concept. In section V we digress from the main theme of this paper and turn 
our attention to the generators of spatial diffeomorphisms of E . Specifically, we show that 
such generators generically do not have a well-defined action on the Fock vacuum. Section 
VI contains a discussion of our results and of open issues. Some technicalities are collected 
in the Appendices. 



6 Here Ej and E^- are slices of constant t, i.e. i = ij and t = tt. Hence r is the radial coordinate on Ej and 
Ej. Since r is in general different from the coordinate R, r = r.i(R) and r = rf(R) express the functional 
relation between the two coordinates on E j and on Sj. 
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In what follows we set H = c = 1. Note that in our discussion of the relation between 
cylindrical waves and PFT, Newton's constant "per unit Z" |4| has been set to ~. 

Note on notations 

The quantum counterpart of a classical quantity, X, is denoted by X, its complex con- 
jugate is denoted by X*. d t X, d r X are denoted by X and X', respectively. The adjoint of 
X is X* and its normal-ordered version is : X :. f(y, k) = 0(1/ k n ) means that, for k — > oo, 
3 a positive number C such that \f(y, k)\ < C/k n , with C independent of y. 



II. EVOLUTION AS A BOGOLUBOV TRANSFORMATION 

The axisymmetric free scalar field, ip(R, T), satisfies the wave equation on the flat space- 
time of equation (|1.2|) . Hence, we have that 

** + sg + m =0 . (21) 



dT 2 dR 2 RdR 
The mode expansion of <fi which is regular at the axis is 

7p(R, T) = — / dk J (kR) [ a(k)e- tkT + a*(k)e lkT } . (2.2) 
v2 Jo 

Consider an arbitrary axisymmetric foliation of the flat spacetime (jl.2j) coordinatized by 
(i, r) (each leaf of the foliation is at constant t; r is the radial coordinate on each leaf). Let 
S be some slice at time t. The momentum, 7r(r, t), (canonically conjugate to the scalar field 
ip{r,t) on this slice) obtained through a Hamiltonian decomposition of the free scalar field 
action with respect to the foliation [3| is 

* = ( fR'-RT') [{R ' 2 " T ' 2) ^ " ^ ~ fT ' W] ■ (2 ' 3) 

Here (R = R(t,r),T = T(t,r)) denote the spacetime coordinates (R,T) evaluated at the 
point (r, t) on the foliation. As is conventional in the Hamiltonian description, we shall 
suppress the label t and denote the scalar field and its momentum on X by (V ; ( r ) ) 7r ( r )) an d 
the restriction of the spacetime coordinates R,T to S by R(r),T(r). A useful form of ir(r) 
for our subsequent considerations is obtained by re-expressing the r, t derivatives of ip in Eq. 
(12. 3|) in terms of R, T derivatives so that 

7r(r)=R(R'^, T +T'^, R ) , (2.4) 

where ip,T and ip,R denote derivatives of ip with respect to T and R respectively. 

Next, we trade (ip(r), 7r(r)) for suitable linear functionals thereof. We require that these 
functionals reduce to the mode coefficients a(k) of Eq. (J2.2)) when S is chosen to be the flat 
slice S at T = with radial coordinate r = R. Note that from Eq. (J2.2|) and Eq. (j2.4j) . the 
field and the momenta on E , i[>o and 7Tq, are 



oo 



ipo = -^= / dk J {kr) [a(k)+a*(k) 
v2 Jo 

i 

75 ./o 



tt = -r / dk k J {kr) [ a{k) - a*{k) ] . (2.5) 
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Using the orthogonality of Bessel functions we can invert Eq. (|2.5)) to get 

1 f°° 

a(k) = —= / dr [i J (kr)ir (r) — r J\(kr) ip (r) ] . (2.6) 
v2 Jo 

The evolution from S to £ of a{k) = a(k)[ipo, vr ] (where the right hand side of the equality 
represents the functional of ip , n Q appearing on the right hand side of Eq. (j2.6j) ) is obtained 
by replacing the fields ipo^o by their evolved images tp(r),n(r) on S. Thus we have that 
(in obvious notation) 



1 f°° 

a s(^) = -j= / dr [i J (&r)7r(r) — r Ji(kr) ip'{r) 
v2 Jo 



(2.7) 



Using Eq. (j2.2|) and Eq. ()2.4|) . ip(r), 7r(r) can be re-expressed in terms of a(k), a*{k) to obtain 
an equation of the form 



dk 2 [a(ki, k 2 ) a(k 2 ) + (3 (hi, k 2 ) a*(k 2 )] 



(2.8) 



Although we do not display the calculations here, it can be checked that a and f3 satisfy the 
standard Bogolubov conditions (see for example 0]). On general grounds this is expected 
since, as mentioned earlier, evolution is a canonical transformation. (3(k\, k 2 ) is given by 



P(h,k 2 



2 



dr r 



Jiih^J^hR) - (-)J (k ir )J (k 2 R) 



R' 



'-)Jo{kir)J x {k 2 R) + J^k^JoihR) 
r 



r 



Jk 2 T 



(2.9) 



Since we do not need the expression for a in this work, we do not display it here. 

From the work of Shale [13|, the canonical transformation defined by Eq. ()2.8|) is unitarily 
implemented iff j3 satisfies the Hilbert- Schmidt condition 



dk\ 



dk 2 |/3| 2 < oo 



(2.10) 



We end this section with a specification of the class of embeddings (T(r), R{r)) for which 
we analyse (J2.1U)) . We require that the embeddings be 'Minkowskian' near the axis of 
symmetry, r = 0, and at spatial infinity. Specifically, we require that R(r) — r and T'{r) be 
compactly supported away from r = and r — > oo. Hence, on every slice £ we require that 
there exist positive a,b,c,d with b > a, d > c such that 



T'(r) =0 for r G [0, a) U (b, oo) 



(2.11) 



R(r)-r = forr€[0,c) U (d, oo) . (2.12) 

Note that (J2.11j) and (J2.12j) ensure that the slice E does not have a 'cusp' at the axis. 
For future reference, we note that the spatial nature of £ implies that 

R' 2 -T' 2 > . (2.13) 

In particular, the above equation in conjunction with ()2.12|) implies that 

R' > , (2.14) 

so that R can indeed be used as a radial coordinate on E. Finally, we shall also assume that 
R, T are C°°-functions of r. 
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III. THE 'HALF PARAMETRIZED' EVOLUTION 



In this section we show that the 'half parametrized' evolution satisfies the Hilbert-Schmidt 
condition, Eq. (j2.10j) . As mentioned in the introduction, 'half parametrized' evolution is 
defined as evolution in which the radial coordinate 'r' on £ is chosen to be the restriction 
of the spacetime coordinate 'i?' to S. Setting R = r in Eq. (j2.9j) . we get 

k f°° 

P{kiM) = -rf / drrdJtih^^k^-Joih^Joifar)} 
z Jo 

~ tiHk^J^hr) + J^hr) J (M] T')e lk2T . (3.1) 

For later calculations, however, a more symmetric form of (3 turns out to be convenient. As 
shown in Appendix IH| using Bessel function identities, -£ [z n J n (z)} = z n J n ^i(z), JjJ (z) = 
—Ji(z), their orthogonality relations, and integration by parts in conjunction with Eq. (|2.11|) . 

we get 

P(k 1 ,k 2 ) = - t - klh I dr riMhr^ihr) + Mk^Mhr)] T'e^ T . (3.2) 
2 [h + k 2 ) J a 

Clearly, divergent contributions to Eq. (j2.10j) may be obtained only when ki or k 2 or both 
diverge. We estimate these 'ultra-violet' contributions to the integral in (j2.10j) by considering 
the following four exhaustive cases. 

A.l. k 2 oo, < fci < k\ /3 . 

A. 2. h -> oo, < k 2 < k\ /3 . 

B. l. k 2 , ki — > oo, k)j 3 < k\ < k 2 . 
B.2. ki, k 2 — > oo, k\^ 3 < k 2 < k\. 



[Case A.l] k 2 -> oo; < h < k 2 ; k 2 a > 1. 

Since k 2 a 3> 1, we use the asymptotic expansion of Bessel functions, Appendix El for the 
ones with k 2 r as their argument. The integral contributing to the first term in Eq. (J3.2j) is 

b 

dr rJ (fcir)J 1 (fc 2 r) T'e ik2T 

= —L= t dr v ^T'J (A; 1 r){e^ / V fc2(r+T) + ^I^-^-t) +0 (1_\\ .(3.3) 
\Z2irk 2 J a \hj 

Since T' has compact support away from the origin, ^T'Jo(kir) is a smooth, bounded 
function of compact support. Further, for space-like embeddings, (r ± T)' 7^ because 
\T'\ < 1 (this follows from Eq. ()2.13j) with R = r). Then the integral in Eq. (j3.3)l is a linear 
combination of integrals of the form considered in Lemma 1 with f 3 = 0. A similar analysis 
yields the same conclusion for the second term in Eq. (j3.2|) . Hence, using Lemma 1, (3 can 
be estimated as 

1^)1 <J^k) 



1.3/2 



Mi + 



.3/2 



Mo 



(3.4) 
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for suitable positive real numbers M±, M 2 - It follows that 

\P{Kk2)\ 2 <%[Ci + C 2 h + CzkJ\ , (3.5) 

K 2 

for suitable Cx, C 2 , C3. This gives the estimate 

,1/3 

jf" *,)!'< l^ + ^-ij + f-ij , (3.6) 

which in turn implies convergence of the integral in Eq. (|2.10j) in the domain of integration 
considered here. 



[Case A. 2] kx — > 00; < k 2 < ki' ; feia > 1. 

The only non-symmetric part (with respect to kx and fo) of /3(kx,k 2 ) in Eq. ()3.2|) is e tk2T , 
which results in the terms e *( /C2T±fclT ') rather than e tk2 ^' ±T ^ in Eq. (j3.3|) . However, Lemma 1 
still applies. The rest of the argument is similar to the one for [Case A.l], and also yields a 
convergent contribution to Eq. (J2.1(J|) in this domain of integration. 



1 /3 

[Case B.l] k 2 , kx — > 00; k 2 < k\ < k 2 ; kxa, k 2 a 3> 1. 

In this case we can expand all the Bessel functions using their asymptotic forms for large 
arguments. Then we have that, 



(3{kx,k 2 



% (kxk 2 y/ 2 

8tt (kx + k 2 ) 



dr 



8 cos(fci + k 2 )r — [sin(A;i + k 2 )r 



+ 2cos(kx-h)r) (- + - )- + 0[- 2 



kx k 2 



IJkiT 



T'e 



(3 + + (3. + K 



(3.7) 
(3.8) 



Here j3 + is the sum of the first two terms (i.e. the cos(fci + k 2 )r and the sin(&i + k 2 )r terms), 
while (3- is the third term (with the cos(fci — kijr). K denotes the remainder of order l/k\. 

First we show that (3+ falls off faster than any power of ^ kl \ k2 ^ ■ The two terms in (3 + are 
of the form 



r b 

f±(h,k 2 ) / drg±(r){e i{kl+k2)r ±e- l{kl+k2)r }T't 

J a 



ik 2 T 



(3.9) 



where g±(r),T(r) are smooth bounded functions of r in [a, b}. Consider the first term in 
()3.9|) . Define x( r ) '■= T(r) + r, t := fci ^ fc „ and Xt as 



ki+k2 

Xt ■■= tr + (l- t)x(r) . 



(3.10) 

Since \T'\ < 1 for a space-like embedding, x't m the domain of the integration. This 
implies that there exists a C°° function ipt(y) = r where Xt{ r ) — V- Then, 



r 

drg±(r)e i( - kl+k2)r T'e ik2T = / 

J Y 



Xt(b) 
Xt(a) 



dyg±{r{y)) d ^T' {r{y)yy^ +k * 
dy 



(3.11) 
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This is an integral of the form in Lemma 2, L.2.a. Therefore the integral falls off faster than 



any power of 



i 



(ki+k 2 ) 



in the large k\ and k 2 limit. Similar arguments apply to the second 



term in (|3.9jl . so we conclude that (3+ falls off faster than any power of ( fcl+fca ) 
hi, k 2 limit. 

Next, consider K. Since k 2 + k\ > k 2 and T' is bounded 



in the large 



\K\ 



?7T (fci + k 2 ) 



dr T'e ik * T A 



< 



C 4 



1/2,3/2 



h 1 h 
h, 2 «,j 



(3.12) 



with a suitable positive C4. 
Finally, consider /3_, 



4 7 r(fc 1 A; 2 ) 1 /2 



c?r cos((/ci — fe)?" 



ife 2 T 



Since the integrand is bounded, a rough estimate for /3_ yields 
Eq.(EHD 



(3.13) 



0((M;2)~ 1/2 ). From 



|/?r < 1/3+1 + 1/3-1 + 1^1 + 2|/3+/3-| + 2|/3 +J ff| + 2|/3_K| 



(3.14) 



Since falls off faster than any power of ( fcl ^ fca ) > so do |/?+| 2 , |/3 + /3_|, and Hence 
these terms give convergent contributions to the integral in Eq. ()2.10j) . 

Next, from Eq. (|3.12j) . and /3_ = 0((k\k 2 )~ 1 / 2 ), we estimate the contributions from \K\ 2 
and 2\/3-K\ as 

dh \K\ 2 < ° 4 



1/3 
2 



1.5/3 
h 2 



and 



pk2 

JkV 3 



k 



4/3 



(3.15) 



(3.16) 



for some appropriate positive constant C5, which lead to convergent contributions to the 
integral in Eq. (j2.10j) . 

So far all the contributions to Eq. (j2.10j) except the one from |/3_| 2 have been shown to 
be convergent. The rough estimate of (3- = O {{kik 2 )~ 1 ^ 2 ^ , used above, however, seems to 
indicate that this contribution diverges. But, as we shall show now, a better estimate of 
yields a convergent contribution. 

We introduce an auxiliary variable a := k K kl . Clearly, < a < 1. /3_ is then given by 



^{hhy/ 2 J a dr r [ 



fe 2 

r 



,ik 2 {T-ar) _|_ e ik 2 {T+ar) 



'] 



(3.17) 



Since \T'\ < 1 and [a, b] is compact, there exists a number m such that \T'\ < m < 1. 

m > fh. Consider Eq. ()3.17|) when a > m. If we naively 



Furthermore, let m := m + 



1 — m 
2 ' 



„2/3x 



estimate the integral in Eq. (j3.17J) as of 0(1), the prefactor, (k^)^ 1 ^ 2 = 0(l/k 2 ~), gives 
divergent contribution to Eq. (j2.10j) . However, since T 1 ± a 7^ 0, the integral falls off faster 
than any power of l/k 2 by Lemma 2 L.2.a. Thus, both terms in Eq. (j3.17|) give convergent 
contributions. 
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The remaining contribution to be estimated is when a < m (That is, (1 — m)k 2 < hi). 
We obtain, for an appropriate constant C, C > 0, for < a < m, 



C 

2 

2 J a J a 



dr x dr 2 T ( ri ) T '( r2 ) e iHnn)-T{r 2 )] [ cosaA . 2 ( ri + r2 ) + cosa^^ - r 2 )] . 



rir 2 



\P-Y < F 

(3-18) 

Since the prefactor is 0(k 2 2 ), we still need the tiniest power of l/k 2 from the integral 
involving cos ak 2 (r\ ±r 2 ) for convergence. The key point is to change variables from (fei, /C2) 
to (a, /C2) and integrate |/3_| 2 over a. The integral of cosk 2 a(ri + r 2 ) over a gives a factor 
of l/k 2 since (ri + r 2 ) > a. Since dk\dk 2 = k 2 dk 2 da, this term contributes J k 2 dk 2 (k 2 3 ) 
which is convergent. 

The only remaining non-trivial contribution is from the term with cosak 2 (ri —r%). From 
Eq. flH.18j) this contribution can be estimated as 



C 

K 2 JO 



da 



< 



r 1 r 2 



C_ 

h 2 J a J a 



dndr2 lMl^M e iUT(n)-T(ru)] 



c 



— T2 I dr + dr^f(r_,r. 



r\r 2 
sin k 2 r-in 



n)-T(ra)] cos ^(n - r 2 ) 

da cos ak 2 (r\ — r 2 



ko.r_ 



T2 I llr - 1 
k 2 



dr_f(r_,r + ) 



r_>0 



sin k 2 r_m 



+T2 / l,r - 

K 2 



k 2 r_ 
sin k 2 r^m 



r-<0 



k 2 r. 



Ii + l 
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where r_ 



r i — r 2, r + := r i + r 2 and /(r_,r_| 



r'(T-i)T'(r 2 ) 
2rir2 



(3.19) 

is a positive function of 



compact support in (ri,r 2 ), and hence in (r_,r + ) plane. Now we estimate //, 



h 



C 
k\ 



dr A 



m 



i/k: 



+ 



l-e 
2 



dr_/(r_,r + ) 



dr_/(r_,r+) 
sin k 2 r_m 



sin k 2 r^m 



k 2 r-in 



+ 



dr_/(r_,r+) 



sin k 2 r_m 



k 2 r. 



k 2 r. 



The last integral doesn't contribute since, for e < 1 and large enough fc 2; /( r -^' ^ ^ 



(3.20) 







because /(r_,r + ) is of compact support. The term 
be bounded because it's maximum is 1. Then, 



sin A:2r_m 



kir—m 



C 



r+<S 



c 

dr—f mSLX -\- -—■> 



dr. 



in the first integral can also 
. i — * 

/max 



2 Jr+<B 



dr. 
i/ fe f r. 



L> E In /c 2 

< 1 

- 1,4 ~ 7,3 



(3.21) 



where / max = max/(r_, r + ), and / = for r + > 5. D, are appropriate constants. By 
redefining r__ to f_ := — r_ we get the same expression for In, and therefore the same 



10 



estimate. This implies that, 

f m In k 

/ da\(3.\ 2 <M— ^ , (3.22) 

JO ^2 

for an appropriate M, which again leads to a convergent contribution. Since all the terms 
in Eq. (j3.14j) are finally shown to give convergent contributions, the integral in Eq. (j2.10j) is 
convergent. 



1 /3 

[Case B.2] k 2 , k\ — > oo; k^ < k2 < k\] kia, k 2 a 3> 1. 

Using methods similar to those in [Case B.l] it can be checked that, once again, the only 
potentially divergent contribution is from |/?_| 2 with /3_ defined as before, by Eq. (j3.13[) . 

To estimate this contribution, we again introduce an auxiliary variable a := ^ > 
(Note that a in the [Case B.2] is different from the a in the [Case B.l] even though we use 
the same symbol.). First, consider the range < a < m (m has already been defined, see 
[Case B.l]). Then, as in the [case B.l] we get, for a suitable constant C, C > 0, 

I/5-I 2 < -S / f dnd ^'^ 1 ^'^ e ik ^ T ^- T ^ )] [cos ak 2 (n + r 2 ) + cos ak 2 (n - r 2 )] . 

k 2 J a J a r l r 2 

(3.23) 

which is exactly the same form as in Eq. (j3.18|) . The rest follows exactly as in [Case B.l] 
and the contribution to ()2.10|) from \(3-\ 2 for < a < m is convergent. 

The only remaining case is when a > m. We define another auxiliary variable a := 
fcl ~ fc2 = !^<x In the range of fci, k 2 considered with a > m, -^-^ < a < 1. This implies that 
0<1 — a < (m + l) -1 - In terms of a we have, 



P_ = \ f l ' dr TL r e «i[(l-«)T(r)+ar] + e i ftl [(l-a)T(r)-ar]l ( 3 . 24 ) 

87T(fcifc2) 1/2 J a r 

where T'/r is a smooth function of compact support. Let f± — (1 — a)T ± ar. ^ 
because |(1 — a)T'| < (1 — a)m < ^-j- < a. Therefore, both terms in the integral fall off 
faster than any power of l/k\ using Lemma 2 L.2.a. Hence, this contribution from /3_ to 
(12.10)1 is convergent. This concludes our analysis of [Case B.2]. 

Since the integral in Eq. ()2.1()j) with f3 given in Eq. ()3.2j) converges in all domains of inte- 
gration, we conclude that /3 for the half-parametrized case is Hilbert-Schmidt. 



IV. ACTION OF SPATIAL DIFFEOMORPHISMS 

As discussed in the Introduction, evolution from So (with radial coordinate 'i?') to an arbi- 
trary axisymmetric slice £ (with radial coordinate V) is a composition of 'half parametrized' 
evolution followed by the action of a spatial diffeomorphism on the final slice. Since 'half 
parametrized' evolution is unitary (see section III), evolution from S to E is unitary iff the 
action of spatial diffeomorphisms on E is unitary. As we show in this section, however, the 
action of a generic spatial diffeomorphism on an arbitrary Cauchy slice is not of Hilbert- 
Schmidt type. Therefore, even though 'half parametrized' evolution is of Hilbert-Schmidt 
type, generic evolution between E and E is not. 
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We proceed as follows. First, we show that the expression for fe) for a diffeomor- 
phism on E is identical to that for f3(k\, k 2 ) for a corresponding diffeomorphism on So- Next, 
we show that the explicit expression for [3(ki,k2) for generic diffeomorphisms of Eo is not 
Hilbert- Schmidt . 

Consider a slice E with coordinate system R (see the remark after Eq. (|2.14j) ). Let 'd' be 
a diffeomorphism from E to itself which maps the point labelled by R to the point labelled 
by r(R). Explicitly, a point with coordinate R = Rq is mapped to a point with coordinate 
R = r(Ro). The action of 'd' on Cauchy data (ip, ir) on E is a linear canonical transformation 
which can be written as 



(ip)d(r) = / drg d (r,r)ifj(f) , (4.1) 
Jo 

POO 

Wd(0 = / drh d (r,r)7r(r) . (4.2) 



Since 'd' drags the fields at R to r(R), we have that g d {f,r) = h d {f,r) = 8(r,R(r)). 

Next, consider the initial slice Eo also with coordinate system R. Let 'do' be the diffeo- 
morphism from Eo to itself which maps the point labelled by R to the point labelled by 
r(R). Clearly, the action of 'do' on data (ipo,7i Q ) is 



Md (r) = / drg dQ (r,r)ijj (r) , (4.3) 
Jo 

POO 

(vrok(r) = / df/i do (f,r)7T (f) (4.4) 



with g do (r,r) = h do (r,r) = S(f,R(r)). Thus, for every 'gP on E, there exists 'do' on E (and 
vice versa) with the same (slice independent) integral kernels. We shall now show that, for 
the Bogolubov transformation of interest, /3(ki, A^) for 'do' is the same as /3(ki, k 2 ) for 'd'. 

First, consider the calculation of j3(ki,k 2 ) corresponding to 'do'. We compute P(ki,k 2 ) 
along the lines described in section II. Denote the image of the mode coefficient a(k) under 
'd ' as a do (k). Then, we have that, 

1 f°° 

a d (k) = —J^ dr [i Jo(kr)(n ) do (r) - r J^kr) (ipo)' do (r) ] 

= — = / dr dr [i J (kr)h do (f,r)7r (f) - r J\{kr) ^ () (f, r)^ (f) ] . (4.5) 
V2 Jo Jo 

Here, we used Eqs. (|4.3|) and (|4.4j) in the second line. Finally, we express (^o^o) i n the 
second line in terms of their mode decomposition, 



V> (f) = / dk[a(k)u(r) + a*(k)u*(f)} (4.6) 
Jo 

POO 

7r (f) = / dk[a(k)v(r) +a*(k)v*(r)} , (4.7) 



where u(f) = l/\f2J G {kf) = u*(f), and v(f) = —i/y/2J {kf) = —v*(f) (see Eq.()2.5|l). This 
allow us to read (3 from 

POO 

«d (^i) = / dk 2 [u{ki,k 2 ) a(k 2 ) + I3(ki,k 2 ) a*{k 2 )} . (4.8) 
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The computation of (3 for an arbitrary slice E is an exactly analogous process. The action 
of £ gP on data (ip, ir) is given by (|4.1j) and ()4.2|) . Here, (ip, n) is the image of initial data on 
So under half parametrized evolution. In section III we showed that the Bogolubov trans- 
formation for half parametrized evolution between a(k) and a^k) was Hilbert-Schimidt. 
It remains to evaluate the action of the canonical transformation (J4.1j) and (J4.2j) on these 
as(fc). From section II, a^(k) and are related by 

1 f°° 

a-s(k) = —= I dr[i J (kr)Tr(r) - r J x {kr) ip'{r) } , (4.9) 
v2 Jo 

and 

POO 

t/j(r) = / dk [a s (Jfe)u(f) + a^(k)u*(r)} (4.10) 
Jo 

POO 

7r(f) = / dk[ax(k)v(r)+ax(k)v*(r)) . (4.11) 
Jo 

From Eqs. (j4~Tj) . p~2j) and the image of a E (A;) under 'd' is 

1 Z" 00 

(a s ) d (A;) = jf rfr [i Ukr){ir) d {r) - r J x {kr) &)' d (r) } 



1 



oc 



= —= / dr dr[iJ (kr)hd(f,r)7r(r)—rJi(kr)g d (r,r)ifj(f)] . (4.12) 
y2 Jo Jo 

Substituting Eqs. (f4~TUl) and (|4~TT|) into Eq.(JHEJ), we can read off (3 from 

POO 

M^i) = / dk 2 [a(k h k 2 ) ax(k 2 ) + P(k u k 2 ) aUh)] ■ (4.13) 



Since hd = hd , gd = gd , (3 is identical to the expression obtained when S = S . Therefore, 
without loss of generality, we shall use (3 corresponding to the action of 'g?o' on the initial 
slice S . 

(3(ki, k 2 ) for the initial slice So can De easily computed by setting T = 0, and T" = in 
Eq.flilJ) so that 

k f°° 

(3{k 1} k 2 ) = ^ dR [rJ.ih^J^hR) - RJo(k ir )J (k 2 R)} . (4.14) 
* Jo 

Note that r Ji{k\r) — RJi(k x R) and Jo(k x r) — Jo(k x R) have support only in [c, <i] since R = r 
outside [c,d], Eq. (j2.12j) . This together with Eq. ()C3j) implies that, for a smooth function 
g(R) with compact support in [c, <f|, where c < c < d < d, such that g(-R) = 1 for R G [c, d], 
we can write /3 as 

f3{h,k 2 ) = yJ 5 dR9(R)[rJi(hr)Ji{k 2 R)-RJo{hr)J {k 2 R)} (4.15) 



2 



dRg(R)R [JifafyhfaR) - Joik^Joi^R)} (4.16) 



Ar,fl](fci,fc2)-Ai«,fl]( fc i> fc 2) • (4.17) 
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To show that the integral in Eq. (|2.10|) with f3 given in (j4.15|) to ()4.17|) diverges, it is 
sufficient to show that the integral diverges over any subset of the domain of integration 
because the integrand in Eq. ()2.10|) is always positive. We choose this subset as follows. We 
change variables from (ki,k 2 ) to (a := ki/k2,k 2 ). Let a Q G R be such that a > 0, ao 7^ 1- 
We shall restrict attention to the subdomain of integration defined by «o — 5 < a < «o + 5 
and k 2 (o) < k 2 < 2k 2 (o), ^2(0) —> 00 • Here 5 is chosen such that 5 <C a so that, in the range 
considered, a > and q^I. 

Similar to the previous section, we use the asymptotic expansion of Bessel function for 
large arguments, Appendix B, to estimate (3. We get, 



dRg(R) 



r + R. (3r - R)(R + ar) 1 



15r -9R)(R 2 + a 2 r 2 ) 




8a^/{rR) 3 
[18r + 2R) 



ko 



a: 



k 2 



ik2(R+ar) 



+ c.c. (4.18) 



1 



47m 1 / 2 



dRg(R) 



r-R (3r + R)(R-ar) 1 

1^- 



'rR 




(15r + 9R)(R 2 + a 2 r 2 ^ 



+ 



%a^(rRf 
(18r-2R) 



a 



^FRJ 



k\ 



e ik 2 (R-ar) _|_ c c 



/?m]+a;. R i+c(^ 



R] 



(4.19) 

(4.20) 
(4.21) 



where 'c.c' refers to 'complex conjugate', is ()4.18|) and A r ™ is f)4.19j) . /^[^r] in Eq. ()4.17|) 
can also be written down in a similar form by replacing r with R in ()4.18j) and ()4.19j) . Similar 
to Eq. (j4.21| ) we can also write /3[ RtR \ as a sum of (3? RR y and a term of 0(1/^) . Then, 

P is, 

' 1 



P = ^R] 



+ P[r,R] + AfLifl + Aitifl + ^ 



(4.22) 



Now, we show that Py RR ^ P[ RR ], and ft^ R ] are of rapid decrease in k 2 . Consider fl^ R ] 
first. The integrals in fiym are of the form in Lemma 2, L.2.a, because R' > and a > 
in the domain of interest implies that 



-(fl + tt r) 



dr , 



A similar argument holds for ft^^. The contribution from /3| 
using Lemma 2, L.2.a. Explicitly, A"^™ is, 



'[R,R] 



(4.23) 

can also be estimated 



dRg{R) 



Ana 1 / 2 

Since a ^ 1, we have that 



[a - 1) . 1 
~2aR~ l ~k 2 



1 3a 2 + 2a + 3 1 



16i? 2 



0^ 



;i-a)i?^0 . 



^2 J 



e ik 2 (l-a)R _|_ c c 



(4.24) 



(4.25) 
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Therefore, (3r R ™ is of the form in Lemma 2 L. 2. a in the domain of interest. 




Since P^rr^, P[rr] an d P^r] are all °f the form in Lemma 2 L.2.a, they are functions of 
rapid decrease in k 2 . Then, from (|4.22j) . we get, 




A rough estimate of P^ r R ^ from ()4. 19|) gives P< r Ri = 0(1). This implies that, 




A term of 0(1/ k^) gives a convergent contribution to the integral in Eq. (|2.1Uj) . and it only 
remains to estimate A~#i with more care. 

As in we show (for a large class of embeddings R(r)) that a judicious use of the 
stationary phase approximation to estimate PZri obtains a divergent contribution to (|2.1(J|) . 
Specifically, we demand that the embedding (and our choice of «o) be such that the following 
"genericity" conditions hold. 

(a) There exists a number ^ 1 with the property that R' = a Q at only finitely many 
points rv, 1 = 1, N in [c, d]. 

(b) R" ^ at 77. 

(c) At least one of /j, / = 1, N, fsee (!4.33[) below) is non-zero. 

Condition (a) implies that there exists a fixed finite number of critical points for the exponent 
in (j4.19|) at a = a - Note that the condition for critical points is R' = a. This condition 
defines a map 'A' from r e S to a G (cko — ^ a o + ^) • From this point of view, (b) states that 
the differential of this map is non-degenerate at r = 77 (with A(rj) = «o)- An application 
of the inverse mapping theorem to A implies that the number of critical points iV does not 
change as a varies in a sufficiently small neighborhood of «o (i.e. if S is chosen small enough) 
and that the location of the critical points varies continuously with a in this neighborhood. 

After these preliminaries, we proceed with our estimate for flZm- The leading order term 
in flUEl is, 



We apply the stationary phase approximation to estimate (|4.28j) . As discussed above, there 
are iV critical points rj where R' = a and we estimate (|4.28j) as 




(4.28) 




(4.29) 



where 




(4.30) 



Gi := arj — R(ri) 



(4.31) 
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Note that 77 depend on a, i.e. 77 = 77(a). 

The stationary phase approximation applied to the subleading order terms in ()4.19|) yields 
an estimate of 0(l/k^ 2 ), which is the same order as the "error term" in Eq. (j4.29|) . Then, 
we have 



a; 



1 



r,R] 47ra i/2 I k 



2n 



1/2 2N 



1=1 



ie 



+ 



k: 



3/2 



(4.32) 



where, for I = 1,2,...,N, fi = fj, G/ = g j, //+jv = //, 7+iV = -Gj l with fi,...,f N and 
G\, ...,Gn defined in Eqs. fl4.30j) and ()4.31j) . It is possible that Gi = Gj for some I ^ J, 

I, J = 1, 2, 2N. By suitably redefining the fj, with the redefined fj denoted by //, the 

sum J^/=i fi£~ lk2Gl can be rewritten as a sum over a finite number N < 2N of terms, each 

of the form fie~ ik2dl , I = 1,2, N, G T ^ Gj for I ^ J. Then, we get, 



Ifi 



- |2 
r,R]\ 



> 



(4.33) 



For generic embeddings, at least one of the fi^O and as mentioned in (c) above, we restrict 
our attention to such embeddings. 

Now, we are ready to estimate the integral in Eq. (j2.10|) in the domain of interest. Using 
Eqs. flOU) and (El tne integral is, 



2/,- 



2(0) 



2(0) 



rao+S 

k 2 dk 2 / da\f3(k 1 , k 2 )\ 2 

J Q.q—8 



1 

> - 



2t 



2(0) 



k 2 dk 2 



2(0) </a -<5 



da|/3i 



ryR]l 



> 



1 



167T 



2t 



2(0) 



2(0) 



ao+S 



ao—S 



da 



a 



16n 



2k. 



2(0) 



1 



2(0) 



2k; 



ao+s da 



A' 



-ifa(Gj-G„) 



„ f ao+5 daC 
dk 2 / — — 
16tt J km J m _ s a k 2 



(4.34) 
(4.35) 
(4.36) 

(4.37) 

(4.38) 



for some appropriate C > 0. Since at least one fj ^ 0, ()4.36|) is equal to AA; 2 (o) for some 
suitably defined positive constant A. (Note that A depends on a and S but is independent 
°f &2(o))- Eq. (|4.37| ) is estimated as of 0(1) (in k 2 (o)) by integrating over k 2 to get 



pceo+5 

L ~5 « 



da ^ _Jlll^_i e -ik2(&i-& 
I^J (Gj - Gj) 



'j) 



2k 



2(0) 



=2(0) 



0(1) 



(4.39) 
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Similarly an integration over k 2 show that ()4.38|) is also of 0(1). All this implies that, for 
large enough fc 2 (o) 

/•2&2(q) r-OL Q +5 a 

/ k 2 dk 2 da\(3{k h k 2 )\ 2 > -k m , (4.40) 

which diverges as k 2 (n\ ~~ > 00 • Thus, we have shown that the contribution to ()2.10j) from 
the domain of interest is (linearly) divergent. Therefore, the action of a generic spatial 
diffeomorphism is not of Hilbert-Schmidt type. 



V. NONEXISTENCE OF INFINITESIMAL GENERATORS OF SPATIAL DIF- 
FEOMORPHISM 

After our analysis of finite spatial diffeomorphisms in the previous section, we now turn 
to a study of infinitesimal diffeomorphisms. Specifically, we are interested in investigating 
if the quantum operators corresponding to generators of finite spatial diffeomorphisms are 
densely defined on the dense set consisting of the Fock vacuum and suitably defined N- 
particle states. Preliminary computations along the lines of suggest that the generator 
of a spatial diffeomorphism is well-defined on this dense domain iff its action on the Fock 
vacuum is well-defined. Hence, we investigate if the action of such a generator is well-defined 
on the Fock vacuum. 

Note that this section constitutes a digression from the main theme of this paper in that 
there is no direct relation between the existence (or lack thereof) of the action of such a 
generator on the Fock vacuum and the lack of unitary implement ability of generic finite spa- 
tial diffeomorphisms. In particular, even if the generator turned out to be a densely-defined 
self-adjoint operator, it would not necessarily imply the existence of a unitary operator cor- 
responding to the finite diffeomorphisms it (putatively) generates because there could be 
operator domain problems in defining ("path-ordered") exponentiation. Conversely, had 
finite diffeomorphisms been unitarily implementable, this would not necessarily imply that 
the infinitesimal generators were defined as self-adjoint operators or, even if they were, that 
the vacuum was in their dense domain. 

Hence, our purpose is merely to initiate an exploration of properties of the generator 
of spatial diffeomorphisms in the hope that such an exploration may be useful for further 
work in the field. The result of our computations is that operators corresponding to generic 
generators of spatial diffeomorphisms 7 of the flat initial slice So do not have a well-defined 
action on the Fock vacuum. The generator of spatial diffeomorphisms of a slice S in classical 
theory is J E drf(r)Trtfj' where f(r) is the radial shift vector. When £ = S , a straightforward 
computation shows that the action of the corresponding normal ordered quantum operator 
on the Fock vacuum state is, 

roo — . lii 

I I til 1 j yr ) . (fyi, u yifj yi, u j . |u/| '"' ' ,? * '* ' '' '' 

JO 

where 



OO POO 

drrf{r) : ^(r, 0)^(r, 0) : |0> f = J / dk x / dk 2 \F{hk 2 )[' , (5.1) 
° Jo Jo 

poo 

F{h,k 2 ) = hk 2 drr/(r)[J (A;ir)Ji(A: 2 r) + Ji(A;ir)Jo(A: 2 r)]. (5.2) 
Jo 

7 The notion of genericity used here is defined later in this section and in Appendix [E] We do not know 
how it relates to the genericity conditions of section IV. 
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Here the radial shift, f(r), is assumed to be of compact support away from r = 0. We 
use Eq. (J2.2J) and the canonical commutation relations between the a(k) and a\k) in the 
left hand side of Eq. (|5.1|) to get Eq. ()5.2|) . If the integral on the right hand side of Eq. ()5.1|) 
diverges, the action of the generator is not well-defined. In the rest of the section we show 
that this is indeed the case. 

To show that the integral in Eq. (J5.1|) diverges, we only need to show that it does so for 
a sub-domain of integration because the integrand is always positive. We restrict attention 
to the sub-domain k 2 > ki, k 2 — k\ — a with a of 0(1) and k\,k 2 — > oo. Using the 
Hankel transformation, G(k), of rf(r) with respect to Jo(kr) (see Eq. (jDl)0 and identities 
in Appendix [Dj we get 



F(ki, k 2 ) = k x k 2 



k 2 +k\ 



dk kG(k) 



k 2 —k\ 



9{k x , k; k 2 ) 9(k 2 , k; h) 



Txko 



irki 



+ h 



k 2 — k\ 



dk kG(k) , (5.3) 



where 9(ki,k;k 2 ) is the angle between k\ and k in the triangle formed by ki,k and k 2 , 
and similarly for 9(k 2 ,k;ki) 8 . By adding and subtracting terms, and using 9(k 2 ,k;ki) + 
9(ki, k; k 2 ) + 9(k 2l k±; k) = ir, we get 



rki+k 2 

F{k u k 2 ) = h dk kG(k) 

Jo 



7T 



J k 2 —k\ 

+ / dkkG(k)(k 2 -kx] 

J k 2 —ki 



6(k 2 , k] ki 



7T 



(5.4) 
(5.5) 
(5.6) 



Before we jump into the rather lengthy estimate of each term, it would be helpful to 
state the main results here. It turns out that, for sufficiently generic embeddings, the 
leading contribution to F(ki, k 2 ) comes from ()5.5|) in the domain we are interested in. Using 
the estimate of (J5.5j) we show that the right hand side of Eq. (J5.1|) diverges. Contributions 
from both (J5.4)) and (J5.6j) turn out to be negligible compared to the one from (J5.5|) . 

First, we estimate the integral in Eq. (j5.4j) . It can be written as 



fci 



ki+k 2 



dk kG{k) = Jfci / dk kG{k) - k x 



dk kG{k) 



(5.7) 



fci+fc 2 



The first integral on the right hand side is zero because, from Eq. (|D2J) . it is equal to rfir) | r= o, 
and f{r) does not have support at r = 9 . Moreover, as shown in Appendix A.l of [5J], G(k) 
falls faster than any power of \jk as k — > oo. Since k > [k\ + k 2 ), k — > oo, and the 
contribution from the second integral is negligible. 

Next, consider Eq. ()5.5|) . There are two contributions, namely, from G(k) and from 
9(kx, k 2 ; k). Since G(k) falls off rapidly when k — » oo, and 9(ki,k 2 ;k) is bounded, large 



8 The notation used here is following. 6(k a ,ki;k c ) is defined as the angle between k a ,kb in the triangle 
formed by k a , kb and k c . 

9 The Hankel transformation with respect to Jo (kr) can be thought of as a two dimensional Fourier trans- 
formation. Since f(r) is Schwartz, so is its Fourier transform, G(k). See Appendix A.l of 
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contributions can arise only when k is small. So, we need to estimate 9(ki, k 2 ; k) when k is 
small. 

Since k\, k2 and k form a triangle, we have in the domain of interest that 



Q 

k 2 = (ki - k 2 ) 2 + 4fcifc 2 sin 2 - 



a 2 + Ak\k 2 sin' 



9 



(5.8) 



Since < 9 < tt, k increases when 9 increases, for fixed k\, k 2 . Therefore, k is small when 9 
is small. For small 9 we can expand sin | in Eq. (j5.8)l to get 



k 2 = a 2 + Akik 2 



fl2 /)4 

--- + 0(9") 
4 48 V ; 



(5.9) 



Now, to see how small 9 and k should be to have a non-trivial contribution, we divide 
the domain into two regions; 9 > (hk^- 1 ' 2 ^ and 9 < (kik 2 y l/2+t where < e < 1/2. 
If = (k^)- 1 / 2 ^ , then k 2 = (hkj 2 * + 0(1). Thus, k ->• oo when k u k 2 oo. Since k 
increases with 9, the contribution from 9 > (kik 2 )~ 1 ^ 2+e is negligible due to the nice fall off 
behavior of G(k) at large k. 

Consider the contribution to Eq. (j5.5|) from 



< k < Va 2 + (M^p 



(5.10) 



This implies that < 4sin 2 #/2 < (kik 2 ) 1+2e so that the upper bound of 9 in this range 
exceeds (hk 2 )- 1/2+e by 0((AaA; 2 )- 3 / 2 + 3e ). Since the contribution from the entire range 9 > 
(k\k 2 )~ l l 2+e is negligible, it makes no difference whether we estimate the contribution from 
9 < (kik 2 Y 1/2+e or from Ea. (l5~TU|) and we use the latter. 

In this range, for large enough ki, k 2 , we have from Eq. (j5.9|) that 



9 



or 



2\ 1/2 



k\k 2 



O 



Then, Eq. (|5.5|) satisfies the bounds 



fc 2 — k\ 



TT 



> 



(kik 2 ) 3 / 2 ~ 3( 



(5.11) 



dk kG(k) 



9(k 2 ,k 1 ;k) 



TT 



ki+k 2 



> 



dk kG{k) 



9(k 2 ,k 1 ;k) 



7[ 
TT 

h 



dk kG(k) 



TT 
k 2 



a 



2\ 1/2 



k x k 2 



^/a 2 + (fcifc 2 ) 2 



dk kG{k)0 



ki+k 2 

dk kG{k) 

h\-\h\ , 



(Ms) 3 / 2 -^ 

9(k 2 ,k 1 ;k) 



(5.12) 
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where we have used (15. lip . As discussed above, is negligible. On the other hand, we 
have that 



y/a 2 +(kik 2 ) 2t 



< 



< 



7r(A;ifc 2 ) 3 / 2 
nih^-se J a 



dk kG{k)0 [(hk 2 ) 3t ] 



= h C 
for appropriate constant C\ and C 2 . Finally \1 2 \ is, 

h 



dkk\G(k)\ 

00 

dkk\G{k)\ 



(5.13) 



7T 



x /a 2 + (feifc 2 ) 2 



1 jkx 

TT\k 2 



1 fki 



dk kG(k) 



( k 2 -a 2 \ 1/2 



cifc kG(k)\[W 



n V k 



|C(a, fci,fc 2 )| 



(5.14) 



If C(a, fei, k 2 ) 7^ as ki, k 2 — > 00 it follows that 7 2 dominates over 7 3 in this limit. The 
leading contribution to ()5.5|) is then ()5.14|) . We now show that, for sufficiently generic shift 
vector f(r) (see Appendix|E]for the definition of genericity), this is indeed the case. The idea 
is to estimate C(a, k\, k 2 ) as a — > and ki, k 2 — > 00. Consider the behavior of C(a, k\, k 2 ) 

< 5 <C 1 and k\,k 2 — > 00. Define <r as a := ak b 2 . Then, 



in the region where a G 



1 2 

/t 2 ^2 



C(a, fci, /c 2 



dk kG(k)Jk 2 - (-^) 2 



dfc kG(k)Jk 2 -(—V 

{kik 2 y 



1,8 ■ 
"2 



(5.15) 

The second integral in Eq. (j5.15|) is negligible since G(k) falls faster than any power of 1/A;[5j. 
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The first integral in Eq. (j5.15j) can be estimated as, 



I 

J — r 



a 



dk kG(k)Jk 2 - (j^) 2 
k 2 



a 



dkkG(k)Jk 2 -(-j) 2 + 

V k 2 



(kik 2 y 



dk kG{k)Jk 2 



a 



01 



0( 



J/2- 



•J/2- 



01 



k. 



8/2 



+ 



+ 



(kik 2 y 



(kik 2 y 



2 

(kik 2 y 



/ik\k 2 y | 
dk k 2 G(k)0(-j) 
k% 
2 



dk k 2 G(k) + 0(-j' 



dkk 2 G{k) . 



(5.16) 



In the second line, we divide the domain of integration because we want to expand 
k 2 — (fx) 2 - This is possible only when k ^> u/k^, i.e. in the second term. Now, 



(kik 2 y 



dk k 2 G(k) = / dk k 2 G(k) 



dkk 2 G(k)- / dkk 2 G(k) . (5.17) 



(kik 2 y 



The second integral in Eq. (|5.17|) is of 0(-fm), while the third integral can be neglected since 

rl 

G{k) falls faster than any power of 1/k, p. This implies that, for a suitable constant E, 



|C(a, A;!, A; 2 )| > 



1 

> 2 



dk k 2 G{k) 



dk k 2 G{k) 



Ea 

1JJ 2 



(5.18) 



provided J °° dk k 2 G{k) ^ 0. In Appendix |E[ we show that this is indeed the case for generic 
G{k). Hence, we conclude that the estimate of Eq. (|5.5|) is, 



k 2 —ki 



TV 



> 



> 



1 fk[ 



2tt V k 2 
1 fk[ 
47T V ko 



dk k 2 G{k) 



dk k 2 G{k) 



D 



(5.19) 



for sufficiently large ki,k 2 . 

Before we estimate E(k 1 ,k 2 ), we need show that Eq. (|5.6|) is negligible compared to 
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Eq. (|5.5|) . For the range of variables we consider, Eq. ([5.6j) satisfies, 



k2+ki 



dk kG{k){k 2 - h 



6{k 2 , k; ki 



k 2 —ki 



71 



rki+k 2 

< a dk \kG{k)\ 

J a 

poo 

< a dk \kG(k)\ 

Jo 

h 2 



(5.20) 



where we used the fact < 9 < tt in the first line. 

Now, we are ready to estimate F(ki, k 2 ). From the Eqs. (|5.7|) . (|5.13|) . f|5. 14|) and (j5.18|) . 



1^1,^)1 > 



> 




dk k 2 G{k) 



dk k 2 G(k) 



(5.21) 



1 2 
77 ? 77 

K 2 «<2 



for sufficiently large k%, k 2 and k 2 — k\ e 

Finally, we have all the necessary ingredients to estimate the integral in Eq. ()5.1|) in the 
range of interest. First, we change variables from (k\, k 2 ) to (a, k 2 ), where a = k 2 — ki as be- 



fore. Then, we integrate |F| 2 over the domain a e 
oo. The integral in Eq. (j5.1j) is, then, 



1 2 
77 5 77 
"a 2 , 



and k 2 G [(fc 2 )o, 2(fc 2 ) ] ] (k 2 )o 



i poo poo 1 

- / dh dk 2 \F(h,k 2 )\ 2 > - 

1 JO JO 1 J(k 2 )o Jl/k? 



2(fe 2 ) r 2/k. 

dk 2 



> 



dk 2 



da \F(ki, k 2 ) 
da{C ) 2 



HCo) 



21-5 



(k 2 )l s [2^ - I 1 "*] - oo 



(5.22) 



for small 5. Therefore, the action of the normal-ordered generator of generic spatial diffeo- 
morphisms of Eq on the Fock vacuum is not well-defined. 



VI. DISCUSSION 

The work of Kuchaf maps the dynamics of cylindrical waves to that of an axisymmetric 
free scalar field along arbitrary axisymmetric foliations of the fixed (2+l)-dimensional flat 
spacetime ([1.2)1 . In this work we studied quantum evolution of this free field operator from 
the initial flat T = slice, So, with radial coordinate R, to an arbitrary slice S with 
radial coordinate r. We showed that operator evolution is unitarily implemented in the 
standard Fock representation when r is chosen to coincide with R on S. The transition to a 
general radial coordinate r(R) from R is obtained by the action of a corresponding spatial 
diffeomorphism on the field variables on E. We showed that for generic choices of r(R), this 
diffeomorphism is not unitarily implemented in the standard Fock space representation. 
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In (l+l)-dimensional parametrized field theory (PFT) where such operator evolution is 
unitary, a Tomonaga-Schwinger type of functional evolution can be defined as the unitary 
evolution of Fock states from the initial T = surface to an arbitrary one [lO] . Our results in 
this work indicate that a similar notion of functional evolution for quantum cylindrical waves 
is not a viable concept for generic choices of foliations starting from the slice S . One may 
inquire if such functional evolution can be defined in the standard Fock representation by 
choosing some other slice than S as the initial slice (this would entail expressing the stan- 
dard mode coefficients a(k) as functionals of data on the new initial slice and analysing the 
Bogolubov transformation corresponding to evolution of these mode coefficients). 10 Though 
we have not investigated this question, we suspect that no such choice of initial slice can 
render generic evolution of operators as unitary. 

Our results do indicate, however, that it should be possible to define functional evolution 
of quantum states along the restricted class of foliations wherein the radial coordinate on 
each slice is R. Such evolution is formally described by the 'half parametrized' formalism of 
It would be of interest to construct the Schrodinger picture states as the unitary image 
of the Heisenberg-picture Fock states and to show that they satisfy a functional Schrodinger 
equation along the lines of [10(, for this restricted evolution. 

Our results have been obtained for foliations satisfying the boundary conditions ()2.11|) 
and ()2.12|) . The fact that T'(r) and R(r) — r vanish at the axis and at spatial infinity 
are direct consequences of the boundary conditions for cylindrical waves in their ADM 
description (see the appendix of |l6|). Our conditions of compact support for these quantities 
are more restrictive than the conditions in It would be of interest to seek, either a 



generalization of our proofs to those boundary conditions, or an alternative set of boundary 
conditions which allows the definition of a consistent Hamiltonian framework and which 
induce conditions on T(r),R(r) such that our proofs still go through. We suspect that 
requiring fall offs of various phase space variables to be faster than any power of r" 1 at 
spatial infinity and any power of r at the axis should suffice for the latter. It would be good 
to check this. 

In section V we showed that a generic infinitesimal generator of spatial diffeomorphisms 
on E does not have a well defined action on the Fock vacuum. What about generators of 
infinitesimal 'half parametrized' evolution on S ? From [^] the relevant Hamiltonian density 
on any slice S with radial coordinate R is given by 

U = \<X - T^y'KR-h - R^T R ^ R f + iityjj . (6.1) 

On £ > T R = 0, it = Rip^T and H equals the standard flat spacetime energy density. 
H(f) := f dRf(R)H is a generator of 'half parametrized' evolution on S . In it was 
shown that for smooth / vanishing fast enough at the axis and at spatial infinity, the 
normal-ordered operator H(f) is densely defined and has a well-defined action on the Fock 
vacuum. Since we know that the standard flat space Hamiltonian, H , is also a well-defined 
operator, we conclude from [Hf that H(f + 1) = H(f) + Ho which generates nontrivial 
evolution at infinity and along the axis is also well-defined. 



Note that an affirmative answer for the existence of generic evolution from the new initial slice would be 
consistent with our considerations here, only if the specific evolution between the new initial slice and Eo 
was non-unitary. 
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It can be checked that the Poisson bracket between two such generators of 'half 
parametrized' evolution on So, H(fi) and H(f2) is the generator of a spatial diffeomor- 
phism. It seems puzzling that despite H(fi),H(f2) being densely defined operators, the 
generators of generic spatial diffeomorphisms do not have a well-defined action on the Fock 
vacuum. A possible resolution is that the relevant shift vector does not satisfy our require- 
ments of genericity since it is built out of f%, f 2 in a specific way. However, we suspect that 
the resolution lies in the well-known problem of domains for unbounded operators (see for 
e.g. [II]): although the vacuum is in the domain of both H(fi) and H(f 2 ), it need not be 
so for their product and/or commutator. It would be of interest to confirm this. 

Next, we comment on the relevance of the recent work of Corichi, Cortez, and Mena 
Marugan ^ to our considerations here. They considered the issue of unitary implement abil- 



ity of evolution in the context of a quatization of the Gowdy model. The dynamics of the 
Gowdy model is also determined by that of a scalar field on an auxilliary 2+1 dimensional 
spacetime. Pierri 0] showed that the scalar field admits a natural quantization based on 
the choice of a specific foliation of the auxilliary spacetime. Corichi et al. and Torre 



21[ showed that in the context of this fixed foliation and Pierri's quantization, scalar field 
evolution is not unitarily implemented. However, in |18| . a quantization was constructed in 
which the evolution of a time-dependent rescaling of the scalar field is unitary. As far as we 
can tell, this quantization is inequivalent to that of 0]. In view of the results of ll|, the 
following questions are of interest in the context of our work here. Given an arbitrary but 
fixed axisymmetric foliation of 2+1 dimensional flat spacetime, is there a quantization in 
which a suitably redefined scalar field variable evolves unitarily? If so, is this quantization 
equivalent to the standard Fock quantization? 

We conclude with some remarks regarding the relevance of our results to a Dirac quan- 
tization of PFT and of gravity. At first glance, our results seem to indicate that a straight- 
forward Dirac quantization of axisymmetric PFT in 2+1 dimensions (and hence of quantum 
cylindrical waves) which is equivalent to the standard Fock quantization, does not exist. 
11 Indeed it seems that the lack of unitary implementability of spatial diffeomorphisms 
for quantum cylindrical waves has an adverse lesson for Loop Quantum Gravity (see for 
example ji^]) since the entire formalism is based on a unitary representation of spatial dif- 
feomorphisms. However, a more careful Loop Quantum Gravity type treatment of PFT 
merits the construction of a kinematic Hilbert space representation for both the matter 
and the embedding variables rather than formally setting the action of X a (x) to be multi- 
plication by the embedding variable X a (x) and that of the conjugate momentum variable 
to be P a (x) = —i sx i( x ) ■ Such a treatment is currently in progress [23j] and suggests that 
when the embedding sector is properly defined in the quantum theory along with an appro- 
priate notion of embedding dependent Hilbert spaces for the matter sector, a satisfactory 
Dirac quantization which is equivalent to the standard Fock quantization can be constructed 
notwithstanding the results of this work and [ill ]. 



11 Note that a classical Hamilton- Jacobi type of canonical transformation along the lines of |9| should result 
in the classical Heisenberg picture constraints which directly yield, upon Dirac quantization, the standard 
Fock space. However the PFT constraints which are naturally available take the form appropriate to 
the classical Schrodinger picture 

HQ. 

Our results imply that a direct Dirac quantization of the latter, 
without further canonical transformations, is not equivalent to the standard Fock representation. 
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APPENDIX A: LEMMAS 



Lemma 1 Consider an integral 

I(p, q) - 



drg(r)fi(pr)e 



iqf2(r) e ipf 3 (r) 



(Al) 



that satisfies 

• p, q > 0. 

• g(r),f 1 (r),f 2 (r),f 3 (r) are C°° . 

• g(r) has support only in [a,bj. 

• ft^Oin [a, b\. 

• fi and all its derivatives are bounded in (0, oo). 
Then, 

\I\ < -Mi + -M 2 , (A2) 
q q 

for suitable constants Mi and M 2 . 

Proof : Let f 2 (r) = V- Since f 2 ^ 0, and / is C°° in [a,b], f 2 has a C^-inverse in [a,b], 
X(y) = r. Then, 



i(p, q) 



dy 



(dx\ 

\dy) 



9 [x(v)} fi \pxiv)} e 



iph [x{y)\ Jw 



(A3) 



where x{°) = a an d xip) — b. After integration by parts we get Eq. (|A2|) because all functions 
and their derivatives in the integral are C°° and bounded. 



Lemma 2 • L.2.a : drg(r)e 



ikf(r) 



faster than for any n ask — > oo, if f'(r) ^ 



in [a,b], g(r) and f(r) are C°° , and g(r) has compact support in [a,b\. 

• L.2.b : J Q b drg(r)e lk f^ — > as | for k — * oo, if f'(r) ^ in [a, 6], g(r) and f(r) are 
C°°, but either g (a) or g(b) or both are non-zero. 
Proof : Similar to Lemma ^ 



APPENDIX B: ASYMPTOTIC EXPANSION OF BESSEL FUNCTIONS 



From jiij pi 99 we have that for large \z\ and \arg z\ < ir, 



Mz) 

JAz) 




2 

HZ 

2 

71Z 



cos(z--)+sm(z--) — 



7T , 



9 



cos u — 



, 37T 3lT 3 , 

cos(^ ) — sm{z — —) h cos(z 

4 4 



4 y 128z 2 
3tt, 15 



+ 0( 



4 y 1282 2 



+ 0( 



(Bl) 

:b2) 



25 



APPENDIX C: DERIVATION OF EQ.rtS^l) FROM EQ. rtCT 

Consider an integral 

POO 

I(k u k 2 ,k) := ik / dr rJ (A:ir)J 1 (fc 2 r)TV fcT(r) 



'//•'■•/uU-i'-)-/i(Mj-[e iir(r) -e i * r W] 



rfr rfMx^r) -MM - Mofar) J (V)][e ifeT(r) - e ifcT (°°)] , (CI) 



where we use integration by parts and Bessel function identities, 4- [z n J n (z)] = z n J n _ 1 (z) 
— JAz) in the third line. Then, 



I{ki,k 2 ,k) + I(k 2 ,ki,k) 

POO 

= (k 1 + k 2 ) I drrlJ^k^J^hr)- J (k ir )J (k 2 r)]e ikT ^ 



5(ki,k 2 ) 5(k 1} k 2 



(h + k 2 ] 


poo 

/ dr 




Jo 

POO 


(h + h] 


/ dr 




Jo 


(h + h] 


e ikT(oo) 



ki 



(h + k 2 ) / drrlJ^k^J^hr)- J (hr) J (k 2 r)}e lkT ^ 
Jo 



Here we use an orthogonality relation for Bessel functions, 

8(k,l) 



dr rJ n (kr)J n (lr) 



k 



(C2) 



(C3) 



in the fifth line. Eq. (|3.1j) is then, 

h 



P(h,k 2 ) 



2(h + k 2 ) 

i k\k '"' 



[I{k u k 2 ,k 2 ) + I{k 2 ,kx,k 2 )} 

dr r[J (fcir)-MM + Ji(hr) J Q (k 2 r)} T'e ik2T 



(C4) 



2 {k\ + k 2 ) J a 

which is Eq. (|3.2|) . In the last line, we use the compactness of T', Eq. (|2.11|) . 
APPENDIX D: DERIVATION OF EQ. (J53|) 

In this section we present the derivation of Eq.f (J5.3|) . The Hankel transform G(k) of 
rf{r) with respect to Jo(kr) is [See, |24j p453] 



G{k) 



dr r rf(r)J Q (kr) 
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(Dl) 



rf(r) 

Using Eq.(jn2) into Ea.(Io~2l we get 



dk kG(k)J {kr) 



(D2) 



F{k 1 k 2 ) = hk 2 / dkG(k)k / dr [J (Ar) Jo(M^i(M + J (Ar) JiC^r) J (fc 2 r)] . (D3) 
'o 



From p411 of |2j 

drJ (kr)J (k 1 r)J 1 (k 2 r) 



6(ki, k; k 2 



rrko 



and 



drJ Q (kr)J Q (k 2 r)Ji(kir) 



1 




if | A; — < A; 2 < |A; + k\ 

if k 2 > k + k\ 
otherwise 



irki 



— if Ifc — k 2 \ < k x < \k + k 2 \ 



if k\ > k + k 2 
otherwise . 

Using Eq.flDU) and Eq.(|D5} into En.(jD"3jl we get Eq.(Q. 



1 

ki 




(D4) 



(D5) 



APPENDIX E: PROOF OF / °° dkk 2 G{k) / FOR GENERIC G(k) 

Recall, from Eq. (jDl|) and Eq. ()D2|) . that G(k) is the Hankel transform of G(r) := rf(r) 
with respect to Jo{kr). Using an integral representation of Jo(kr) and rotational symmetry 
of G(r), G(k) := G(k) can be shown to be a two-dimensional Fourier transform of G(f) := 
G(r)[5|. That is, 



G(k) = / dx 2 G(x u x 2 )e iklXl+lk2X2 



(El) 



However we can also consider G(k\, k 2 = 0) as a one-dimensional Fourier transform of g(xi) 
as defined below, 



G(k 1 ,k 2 = 0) 



/2tt 
1 



dxie iklX1 



"Itx J 



-oo 
oo 



2lT 



dx 2 G(xi, x 2 ) 



dx!e iklXl g(xi 



(E2) 



Note that g{x\) exists and posses good fall-off and smoothness property inherited from 
G(xi,x 2 ). Then, G(ki,k 2 = 0) := g{ki) satisfies, 



'2k 



dhe'^gih 



(E3) 
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Now, 



which implies 



dx\ 
d 2 g(x 1 ) 



'2ti J- 



doc-^ 



xi=0 



2,71 



dk x e lk ^k{g{k x ) 



dkiklg(ki) 



(E4) 



d hk\g{ki) 



(E5) 



Here we have used the fact that G(k 1 ,k 2 = 0) := g(ki), and G(k) = G(||A;||) due to the 
rotational symmetry. 



To prove J °° dkk 2 G(k) ^ for generic G(k) we only need to prove that 



d 2 g{xi) 



J. 1 1 



generically. For fixed X\, we change variables from (xi,^) to (x,\,r)\ x 2 i— ► 
Then, from Eq. ()E2|) and rotational symmetry of G(x), 



xi=0 
2 

X2 = r. 



2 rX 

7T ./,, 



dx 2 G(x 1 , x 2 ) 
G(r) 



drr 



.'■1 



yV 2 — x\ 



(E6) 



Since Gir) is a function with compact support away from r = 0, there exists a constant r 
such that G{r) = for r < r . Then, for xi < r , we get, 



That is, for x\ < r , 



2 rX 

7T ., r 



dx\ 



drr 



Gir) 



— xt 



2 ' x 
7r -- 



drr 



G(r) 



(E7) 



— it 



(E8) 



Since we are interested in the limit X\ — > 0, we can use Eq.(jE8J). Then, 

d 2 g(x 1 ) 



doc 



xi=0 



2 x 

71 



dr 



G(r) 



''0 

00 

dr 1 





r 



(E9) 
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where we use the fact f(r) = for r < ro in the last line. For generic embedding the integral 
in the last line does not vanish, therefore, J °° dkk 2 G(k) ^ 0. 
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